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Abstract 

The prospect of detecting relic gravitational waves (RGWs), through their imprint in the cos- 
mic microwave background radiation, provides an excellent opportunity to study the very early 
Universe. In simplest viable theoretical models the RGW background is characterized by two pa- 
rameters, the tensor-to-scalar ratio r and the tensor spectral index Uf. In this paper, we analyze 
the potential joint constraints on these two parameters, r and nj, using the data from the up- 
coming cosmic microwave background radiation experiments. Introducing the notion of the best 
pivot multipole we find that at this pivot multipole the parameters r and are uncorrelated, 
and have the smallest variances. We derive the analytical formulae for the best pivot multipole 
number i^, and the variances of the parameters r and n^. We verify these analytical calculations 
using numerical simulation methods, and find agreement to within 20%. The analytical results 
provides a simple way to estimate the detection ability for the relic gravitational waves by the 
future observations of the cosmic microwave background radiation. 
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I. INTRODUCTION 



lenges for current and future cosmic microwave background radiation (CMB) ex- 
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. The RGWs are produced in the early Universe due to the superadiabatic 



Detection of rehc gravitational waves (RGWs) can be arguably considered one of the most 
important cha 

nr 

periments |2 

amplification of zero point quantum fluctuations of the gravitational field [6]. For this reason, the 
RGWs carry invaluable information about the early history of our Universe that is inaccessible to 
any other medium (see review |7l] for detailed discussion). 

A whole range of scenarios of the early Universe, including the inflationary models, generical 



predict a RGW background with a power-law primordial power spectra 
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In fact, the existence of RGWs is a consequence of quite general assumptions. Essentially, their 
existence relies only on the validity of general relativity and basic principles of quantum field 
theory ^. The RGW backgrounds are conventionally characterized by two parameters, the so- 
called tensor-to-scalar ratio r and the primordial power spectral index of RGWs nt (explained in 
detail below). 

The RGWs leave well understood imprints on the anisotropics in temperature and polarization 

More specifically, RGWs produce a specific pattern of polarization 



of CMB 
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in the CMB known as the i?-mode polarization [l5|. Moreover, RGWs produce a negative cross- 



correlation between the temperature and polarization known as the TE'-correlation 
"3, 



(see also 



181, 



201, 



21 
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24l]). The theoretical analysis of these imprints along with the data from CMB exper- 



iments allows to place constraints on the parameters r and nt describing the RGW background. 

The current CMB experiments are yet to detect a definite signature of RGWs. It is hoped that, 
in the near future, with the launch of the Planck satellite Q] together with a host of ground-based 

and balloon-borne 4| CMB experiments as well as the proposed satellite mission CMBPol Q], 
we shall be able to detect a definite signature of the RGW background. In light of this prospect, it 
is important to be able to effectively constrain the parameters r and nt- A number of papers have 
discussed the current and potential constraint on the tensor-to-scalar ratio r 25[. However, most of 



these works either ignore the constraint on the spectral index nt, or make simplifying assumptions 
about its value. One of the common simplifying assumptions is the so-called "consistency relation" 



nt 



-r/8 |2e 



It should be noted that, the consistency relation is valid 



models of inflation namely the single-field slow-roll inflationary model 



only in the simplest 



23]. For a detailed 
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critical discussion of inflationary predictions and data analysis based on these predictions see [29|. 
In order to keep our discussion sufficiently general we shall not use this consistency relation in our 
analysis. 

The constraints on the parameters r and rit, characterizing the RGW background, will give us 
a direct glimpse into the physical conditions in the early Universe. In particular, they will allow to 
place constraint on the Hubble parameter of the early Universe [sfll] , which in the case of inflationary 



2Q|. More specifically, 



models would correspond to the constraints on the energy scale of inflation 
the amplitude of the RGW power spectrum at a particular wavelength, characterized by r and 
rit, determines the Hubble parameter at the time when the particular wavelength left the horizon. 
Thus, the determination of r and nt would give a direct measurement of the time evolution of the 
early Universe, and provide an observational tool to distinguish between the various inflationary 
type models. In addition, the spectral index nt has a special character if the RGW background 
are generated in a primordial Hagedorn phase of string cosmology [31] or inflation in the loop 
quantum gravity [3^] , so the determination of nt provides an observational way to test or rule out 
these models. 

In this paper we shall analyze the joint constraints on two parameters r and nt that would 
be feasible with the analysis of the data from the upcoming CMB experiments. In general, there 
will be a non- vanishing correlation between parameters r and ut [^,[3^. As will be explained in 
the following sections, the definition of r and nt depend on a reference scale characterized by a 
multipole number Iq, which may be chosen arbitrarily. We shall show that with an appropriate 
choice of this multipole number, which we shall call the best multipole number £j (following the 



terminology of [3j]), the parameters r and nt become uncorrelated and have the smallest possible 
variances. We shall derive approximate analytical expressions for the variances and the correlation 
coefficients, followed by an analytical calculation of the pivot multipole i*. Using the Markov Chain 
Monte Carlo (MCMC) simulation methods, we shall verify our analytical results and evaluate the 
expected constraints for realistic CMB experiments. 

The outline of the paper is as follows. In Section II we shall introduce and explain the notations 
for the power spectra of gravitational waves, density perturbations and various CMB anisotropy 
fields and briefly explain how they are calculated. Furthermore, in this section we shall explicitly 
state the simplifying assumptions that we shall be using throughout the paper, and explain the 
limits of their applicability. Following this, in Section III, we shall calculate analytically the 



expected variances and the correlation associated with the parameters r and rif. We shall show 
existence of the best pivot multipole scale for which the variances of the corresponding r and 
nt are minimal and the correlation between them vanishes. In Section IV we shall confirm our 
analytical results using numerical calculations. Finally, Section V is dedicated to a brief discussion 
and conclusions. 



II. POWER SPECTRA OF COSMOLOGICAL PERTURBATIONS AND CMB 
FIELDS 



The main contribution to the observed temperature and polarization anisotropies of the CMB 
comes from two types of the cosmological perturbations, density perturbations (also known as the 



IC 
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14, hj. These 



scalar perturbations) and RGWs (also known as the tensor perturbations) 
perturbations are generally characterized by their primordial power spectra. These power spectra 
are usually assumed to be power-law, which is a generic prediction of a wide range of scenarios of 
the early Universe, including the inflationary models. In general there might be deviations from 
a power-law, which can be parametrized in terms of the running of the spectral index (see for 



example 



27l|). but we shall not consider this possibility in the current paper. Thus, the power 



k ^ 



spectra of the perturbation fields have the form 

Pn{k) = Asiko) (^^J , (1) 

P,{k) = Mho) ^-J , (2) 

for density perturbations and the RGWs respectively. In the above expression ko is an arbitrarily 
chosen pivot wavenumber, is the primordial power spectral index for density perturbations, 
and nt is the primordial power spectral index for RGWs. ^<j(^o) and At{ko) are normalization 
coefficients determining the absolute value of the primordial power spectra at the pivot wavenumber 
k(). The choices of = 1 and nt = correspond to the scale invariant power spectra for density 
perturbations and gravitational waves respectively. The quantity P-jiik) is the primordial power 
spectrum of the curvature perturbation IZ in the comoving gauge, i.e. P'R.{k) = k'^ {\lZk\'^) /2tt'^ (see 



35[ for a detailed exposition). The quantity Ph{k) is the primordial power spectrum of RGWs 



and gives the mean-square value of the gravitational field perturbations, in a logarithmic interval 



4 



of the wave-number k, at some initial epoch when the wavelenghts of interest are well outside the 
horizon. 

The relative contribution of density perturbations and gravitational waves is described by the 
so-called tensor-to-scalar ratio r defined as follows 

= . ,.y (3) 

Note that, in defining the tensor-to-scalar ratio r, we have not used any inflationary formulae which 
relate r with the physical conditions during inflation and the slow-roll parameters (see for example 



2g]). Thus, our definition depends only on the power spectral amplitudes of density perturba- 
tions and RGWs, and does not assume a particular generating mechanism for these cosmological 
perturbations. 

Assuming that the amplitude of density perturbations As{kQ) is known, taking into account the 
definitions ([2]) and (l3|), the power spectrum of the RGW field may be completely characterized 
by tensor-to-scalar ratio r and the spectral index nt. The RGW amplitude At (ko) = r(fcn)As(fcn) 
provides us with direct information on the Hubble parameter of the very early universe [30| ■ More 
specifically, this amplitude is directly related to the value of the Hubble parameter _H si^ a time 
when wavelengths corresponding to the wavenumber k^ crossed the horizon ^, ^, 
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Mpl JT 



ko/a=H 



where Mpi = l/VSirG is the reduced Planck mass. 

It is important to point out that, for spectral indices different from the invariant case (i.e., when 
Us ^ 1 and rit ^ 0), the definition of the tensor-to-scalar ratio depends on the pivot wavenum- 
ber kQ. If we adopt different pivot wavenumber ki, the tensor-to-scalar ratio at this new pivot 
wavenumber r{ki) is related to original ratio r(A;o) through the following relation (which follows 
from the definitions ([I]), ([2]) and ([3])) 

Let us now turn our attention to CMB. Density perturbations and gravitational waves produce 
temperature and polarization anisotropics in the CMB characterized by the four angular power 
spectra Cf^, Cf^, Cf^ and Cj^ as functions of the multipole number Here Cj'^ is the 
power spectrum of the temperature anisotropics, Cf^ and Cf^ are the power spectra of the so- 
called E and B modes of polarization (note that, density perturbation do not generate i?-mode of 



polarization [3]), and Cj^ is the power spectrum of the temperature-polarization cross correlation. 
In what follows, we shall use the short hand notations Cj, Cf, Cf and Cf to denote these spectra. 
In general, the various power spectra Cj (where Y = T, E, B or C) can be presented in the 



following form 



Y 



t,ti 



(5) 



where Cj^ is the power spectrum due to the density perturbations (scalar perturbations), and CJ^ 
is the power spectrum due to RGWs (tensor perturbations). 

In the case of RGWs, the various CMB power spectra can be presented in the following form 



for Y = T,E,B, 



(6) 



Similar expressions hold in the case CMB anisotropics due to density perturbations with a single 
exception. Density perturbations do not produce the B-uiode of polarization [iS^. Thus, the CMB 
power spectra have the form jiol 



The transfer functions A^"^'*(A;) (see 



, for y = r, E, 



(7) 



18| for details) in the above expressions translate the 
power in the metric fluctuations (density perturbations or gravitational waves) into corresponding 
CMB power spectrum at an angular scale characterized by multipole 1. In general, these transfer 
functions are peaked at values i ~ (1.35 • 10^ Mpc) x k, which is a reflection of the fact that metric 
fluctuations at a particular linear scale k~^ lead to CMB anisotropics predominantly at angular 
scales 9 ^ kD (where D is the distance to the surface of last scattering). In this work, for numerical 
evaluation of the various CMB power spectra due to density perturbations and gravitational waves, 



we use the publicly available CAMB code 



Since we are primarily interested in the parameters of the RGW field, in the analytical and 
numerical analysis below we shall work with a fixed cosmological background model. More specif- 
ically, we shall work in the framework of ACDM model, and keep the background cosmological 
parameters fixed at the values determined by a typical model [38] 



h = 0.732, nbh"^ = 0.02229, l^^/i^ = 0.1277, 17^ = 0, Tr. 



0.089. 



(8) 



Furthermore, for density perturbations, we shall use a model with primordial scalar perturbation 
power spectrum characterized by an amplitude and spectral index 

As = 2.3 X 10"^ Us = 1.0. (9) 

In light of the above, CMB power spectra produced by RGWs depend on the tensor-to-scalar 
ratio r and the spectral index n^. In general, this dependence is complicated and requires nu- 
merical calculations. For analytical calculations in Section IIIH we shall use a simple analytical 



m 



approximation for this dependence (see for example 

£ \ nt-ht 



= Cj, (^) exp Km -ht)in (f/f„)] . (10) 

Here Cj^ = CY^{r = f,nt = fit) are the spectra calculated for values of tensor-to-scalar ratio and 
the spectral index fixed at fiducial values f and n^, and (.q is the pivot multipole. The approximation 
PU]) can be further simplified, for values of spectral index rit sufficiently close to the fiducial value 
ht (such that {rit — ht)\Q.{l/ Iq) <^ 1) 

Clt ^ Cl (J) [1 + {nt - ht) In {e/£o)] . (11) 



The pivot multipole is closely related to the pivot wavenumber /cq- The approximation (jlOp 
can be derived from ^ and ([6]) under the assumption that the wavenumber k and multipole i are 
linearly related, i.e. k/ko ^ i/^o- This assumption is justified due to the peaked nature of transfer 
functions Ay^k) entering ([6]). Numerical evaluations show that the pivot multipole is related to 
pivot wavenumber by 

io^kox lO^Mpc. (12) 

For illustration, in FIG. [1] we plot the power spectra Cj^ for different value of the spectral index 
Tit- The pivot wavenumber is taken to be /cq = 0.05Mpc~^. As expected, in all the panels the 
spectra with different values of rit converge at £ ~ 500, which is consistent with the prediction of 
the relation ([12]). 

The CMB power spectra Cj are theoretical constructions determined by ensemble averages over 
all possible realizations of the underlying random process. However, in real CMB observations, we 
only have access to a single sky, and hence to a single realization. In order to obtain information on 
the power spectra from a single realization, it is required to construct estimators of power spectra. 
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FIG. 1: The CMB power spectra due to RGWs for various values of the spectral index n^: Cj^ 
(left upper panel), Cf^ (right upper panel), Cf^ (left lower panel), Cf^ (right upper panel). The 
pivot wavenumber is chosen /cq = O.OSMpc^^ (in all the panels), and the power spectra are shown 
for multipoles around the value of the corresponding pivot multipole H. ~ 500. 

In order to differentiate the estimators from the actual power spectra, we shall use the notation 
DJ to denote the estimators while retaining the notation Cj to denote the power spectrum. It is 
important to keep in mind that the estimators DJ are constructed from observational data, while 
the power spectra Cj are theoretically predicted quantities. The probability density functions 
(pdfs) for the estimators are described in detail in Appendix |X1 In what follows, we shall require 
the data from all the power spectral estimators, i.e. DJ for Y = T, E, B and C. Let us denote this 
set of estimators (which we shall sometimes refer to as the sample) as 

{Dj} = {Dl\Y = C, T, E,B;£ = 2,3,- ■ •, Cax}- 

To simulate an experiment, we shall randomly draw a data set from the pdf (|Aip . In 
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calculating the pdf ()Aip . along with parameters given in ([8]) and Q, we set the value of the RGW 
parameters as 

r = f, nt = fit. (13) 

We shall refer to f and fit as the parameters of the input model. 

For analytical evaluations in Section HIH we shall work with Gaussian approximation to the exact 
pdfs JXl]) for the estimators [DJ]. The G aussian approximation is characterized by corresponding 
mean values and standard deviations l^I 

{Dj) = CY, iY = T,E,B,C), 

^ -{Cf + NYW^^), {Y = T,E,B), (14) 



a nY 



D} y (2<'+l)/,ky ' 



a nC 



ic^r+icT+Nlw,-^)(cf+Nfw,-'') 



DV \l (2£+l)/,ky 

Note that, the above expressions for mean values and standard deviations follow from the exact 
pdfs considered in Appendix [Al In the above expression, N'[ are the noise power spectra, /sky is 
the cut sky factor, and We is the window function. 

In the case of the Planck mission [l], considering the channel at 143GHz (which has the lower 



21, 



the noise power spectra, the cut sky factor and 
40| for further explanations) 
^ = 5.58 X 10-VK^ 

(15) 
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foreground level and lowest noise power s pec tra 
the window function are given by jl| (see 

A^f = 1.53 X lO^^/xK^ Nf = 

/sky = 0.65, Wi = exp 
where ^fwhm = 7.1' is the full width at half maximum of the Gaussian beam. 

In this paper, along with predictions for Planck, we shall consider an idealized situation with 
no instrumental noise, full sky coverage and an idealized window function Wi = 1. For this case, 
we shall assume that the only source of noise comes from contribution of cosmic lensing to the 
-B-mode of polarization. In this case the noise spectrum for the i?-mode is close to white with a 
value A'f ~ 2 X lO^^^^K^ [41, 42]. A number of works have discussed methods to subtract the 



lensing i?-mode signal (see for example 



41 



431). In 



43t |. the authors claimed that a reduction in 



lensing power by a factor of 40 is possible using approximate iterative maximum-likelihood method. 
For this reason, as a further idealized but feasible scenario, we shall also consider the case with 
reduced cosmic lensing noise ~ 5 x lO^^f^K^. Thus in the two described examples the noises 
are 

Af (lensing) = 2 x IQ-^fiK^, Af (reduced lensing) = 5 x 10"^^K^; 
Nl = Nf = 0; /sky = 1; Wi = 1. 



(16) 
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Note that, in addition to the instrumental noises and lensing noise, various foregrounds, such 
as the the synchrotron and dust, significantly contaminate the CMB signal. However, it is hoped 
that, using multifrequency observations together with ingenious foreground subtraction techniques. 



future experiments would be able approach the ideal limit of expression (16) (see for instant 44i|). 

Before proceeding, let us briefly mention the notational conventions used in this paper. The star 
superscript denotes the quantities evaluated at the best pivot multipole i*. The hat superscript 
indicates the parameters of the fixed (input) cosmological model, that are used to generate the 
simulated observational data. The summation (product) symbols with subscript ^ or y indicate 
summation (product) for i = 2, ...,£max and Y = C,T,E,B respectively. In numerical evaluation 
we set £max = 1000. 



III. ANALYTICAL APPROXIMATION 

In this section, we shall derive analytical expressions for the estimation of the parameters r and 
rif, the associated uncertainties Ar and Ant and the correlation between these parameters. We 
will show the existence and explain the significance of the best pivot multipole i*. Introducing the 
tensor-to-scalar ratio r* defined at the best pivot multipole, we shall show that this parameter can 
be determined with the smallest possible uncertainty, and is not correlated with the spectral index 
rif. Based on this analysis, we shall discuss the signal-to-noise ratio and detection possibilities for 
various CMB experiments. 



A. Approximation for the likelihood function 

In order to estimate the parameters r and nt characterizing the RGW background, we shall use 
an analysis based on the likelihood function 



451, 



46(1 . The likelihood function is just the probability 



density function of the observational data considered as a function of the unknown parameters 
(which are r and nt in our case). Up to a constant, independent of its arguments, the likelihood 
function is given by 

C = llf{D^,Dj,Df,D^), 



where the function f{Df,Dj,Df,D^) is explained in detail in Appendix Rl 
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For analysis in this section, we shall use a Gaussian function to approximate the pdf of the 
individual estimator DJ, and ignore any possible correlation between different estimators. In this 



case the approximate likelihood function can be written as (see 



401] for details) 



The parameters r and rit enter the above expression through the quantities Cj and ctj^y. In our 
analytical considerations we shall make a further simplification. We shall assume that aj^r entering 
p!7|) is weekly dependent on parameters r and rit and assume a^v = Oj^y (for a justification of 



22 



4011). With this assumption, the likelihood function can be rewritten as 



this assumption see 
follows 

^ y \ J 

In the likelihood analysis, we shall assume that the value of the sought parameter rit is sufficiently 
close to the input value ht. In this case, inserting (fTTj) into (fT8|) . using ([5]), we can rewrite the 
likelihood in the form 

) {l + {nt-ht)h) -4 1^ 



21n/: = 5^5;{ar [(J)(l + (n,-n,)M] -df} , (19) 



where 



'^Dj V^O/ (^Dj 

Note that, in the above expression, the dependence on the data (on the estimators DJ) is solely 
contained in the term cUf . Furthermore, of , and are independent of the RGW parameters r 
and nt- The dependence on r and nt takes a particularly simple form and is contained within the 
square brackets on the right side in ([T9|). 

In order to proceed, it is convenient to introduce new variables 

i = r/f C = {nt - nt){r/f), (21) 

in place of r and rit- In terms of these variables, the likelihood (jl9p can be simplified as 

- 2 ln£ = ^ ^ [af (e + C^.) - dj] ' . (22) 
e Y 

Note that, the dependence on the sought for parameters r and nt, in the above expression, is 
contained in the variables ^ and ^. After a straight forward manipulations (|22|) can be rewritten 

11 



as 

-21n£ = e2(E,Ey«P) + C'(E.Ey(arM') + 2eC(E£EyaPM 

-nEi Ey ajdj) - 2C(E^ Ey ^djbi) + Ee Ey 4^. 
This expression can be rewritten as of 

-21n£ = (E,Er -P) (« - + (E^ErK^f^)^) (C - g:g;g>':| )' 

+2«(E,EKor''() + c, 

where C is a constant, independent of r and nf. This constant is responsible for the overall 
normalization of the likelihood function and will not participate in estimation of parameters. In 
the following subsection we shall use the approximation ([23|) for estimating the parameters r and 
nt- 



B. Posterior pdf and the best pivot multipole 



1. Posterior pdf 



arobability density 



45 



46| 



The constraint on the parameters r and nt, are determined by the posterior 
function P{r,nt). This posterior pdf is related to the likelihood function C by 

P{r,nt) = f{r,nt)C, (24) 

where f{r,nt) is the prior probability density function of the parameters r and nj. In this paper, 
we adopt a flat prior, i.e. 

f{r,nt) = l. (25) 

Thus, in this case, the posterior pdf P{r, nt) becomes equal to the likelihood. Using the approxi- 
mation (j23p for the likelihood, we obtain 



2 In P(r,nO = E «^') " S^^) ^ + E 



E<? Ey a 



+2ec(EE"^'''^^)+^- (26) 

e Y 

The parameters r and nt enter the above expression through the variables ^ and C- For this reason 
it is convenient to firstly consider the posterior pdf for variables ^ and ^. It will be seen that the 
posterior pdf for these variables will have a particularly simple form, namely a bivariate normal 

function. The posterior pdf P(i^, is related to P{r,nt) in the following manner 

d{r, nt) 



PitC) 



5(e,c) 



P{r,nt) = -P{r,nt), 



(27) 
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where 



d{r,nt) 



denotes the Jacobian of the transformation between the two sets of variables calcu- 
lable from (j2ip . For simplicity and clarity, let us firstly consider the constraints on the parameters 
and (. Following this, we shall return to the discussion on r and nt using relation (I27jl . 
Introducing the notations 



= {Ei Ey 4dJ)/{Ei Ey 4') > 6 ^ i/JEi Ey « 



Cp^{EiEY4dJh)/{EeEYi'^Jbir), O^l/yjEiEYi'^Ibiy- 
From (j26p and (j27p . we obtain that expression for the posterior pdf of ^ and ( in the form 



(28) 



re 



c 



^(e,C) = ^exp 
4 







\ iC-Cp?] 




[ my \ 


exp 


[ my \ 


exp 



I Y 



(29) 



2. Best pivot multipole i* 

Let us concentrate on the posterior pdf (j29p . As can be seen, there is a non- vanishing correlation 
between the parameters ^ and C in the case when EiEy '^J'^^e 0. From the definition (j20p . it 
follows that the terms b£ depend on the arbitrarily chosen pivot multipole Iq (corresponding to the 
pivot wavenumber kQ through the relation (jl2p ). For this reason, we can select the pivot multipole 
Iq = £* so as to require 



(30) 



Y 



where 6^ = be\^^^g*. With this choice of pivot multipole i^, and the corresponding pivot wavenum- 
ber fcj , the variables ^(A;^) and C{k^) will have no correlation. We shall refer to this pivot multipole 
£i as the Best Pivot Multipole number. From definitions (j20p of aj and 6^, along with expression 
(jl4p . it follows that the precise numerical value of the best pivot multipole number depends on 
the input cosmological model characterized by the ([8|), Q and (jl3p . as well as the specifics of the 
CMB experiment characterized by noise power spectra, cut sky factor and window function. We 
shall discuss this dependence in more detail below. 

Setting the value of the pivot multipole io = if^ so as to satisfy (pOj) . we arrive at a simplified 
form for the posterior pdf 



(fe'^)-^exp 



ic-Cpy' 




{c-c;y' 


msy 


exp 


my 



(31) 



As a reminder let us point out that, in the above expression, as well as in what follows, we have 
used notations r*, (* and 6! to denote the corresponding quantities calculated for the pivot 
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multipole chosen at the best pivot multipole value i*. Note that for the spectral index of RGWs 
we shall retain the notation nt, since it does not depend on the choice of the pivot multipole. 



3. Constraints on parameters ^* and (* 

Equipped with the posterior pdf (j3ip . let us analyze the uncertainties in determining the pa- 
rameters ^* and . For simplicity of analysis we shall assume 

Cp » C- (32) 
In Section IIIIFl we shall show that this constraint corresponds to a condition that the signal-to- 
noise ratio is large, i.e. S/N S> 1. Taking into account this condition, the posterior function 



may be further approximated in the following manner 



re 



*\2 



exp 



(33) 



2te)2 

Note that, in the above expression, the factor in front of the exponent re^ /^* now becomes a 
constant independent of ^* and Thus, the posterior pdf P{i*,C*) becomes a bivariate normal 
(Gaussian) function for variables ^* and C,* . The position of the maximum and the standard 
deviation associated with the posterior pdf P{C ■, C*) are given by 

CML = ep, Ar = c, cXfL = c;, ^c = c- m 

In the above expression, subscript "ML" stands for "maximum-likelihood", since the maximum 
of the posterior pdf coincides with that of the likelihood function due to (|24|) and (j25|) . Following 
the maximum likelihood parameters estimation procedure, we shall accept the values Cajl ^^^^ 
Cml a^ estimators for the corresponding quantities and C*- It is worth mentioning that, for 
the posterior pdfs considered in this work, the maximum-likelihood values coincide with the mean 
values of the corresponding posterior pdfs. It is worth mentioning that, the assumption ^* ^ ^* 
(which was used to derive the pdf ([33]) ) is equivalent to the requirement » A,^*. 

Proceeding further, we can calculate the correlation coefficient for variables ^* and Let us 
firstly define the covariance in the following manner 



cov{x,y) = {x-x){y-y), (35) 
where the overline indicates averaging over the corresponding posterior pdf. The correlation coef- 
ficient can now be explicitly calculated to give 

^ ^°^(^*^C*) = 0, (36) 

y^cov(^*,^*)cov(C*, C*) 
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as expected, the correlation between the variables ^* and vanishes. 

Taking into account ([28j) . we find that the mean values ^* (also Cml) ^^'^ C* (also Ca/l) depend 
on data {DJ} through quantities dj . However, in our approximation, the standard deviations A^* 
and A(* are independent of the data. They depend on the input cosmological model determined 
by ([8]), Q, (|13p . along with noises, cut sky factor and the window function characterizing the CMB 
experiment. 



C. Constraints on parameters r* and rit 

Let us now return to the parameters of our direct interest, namely the tensor-to-scalar ratio 
(determined at the best pivot wavenumber) r* and RGW primordial spectral index n^. These 
parameters are related to ^* and C* through relations 

r* = f*e, nt = nt + C/C, (37) 

which follow from (j2ip . Taking into account the fact that the quantity ^* is peaked at ^^^^ which 
is sufficiently close to the input value and A^*/^|^^ <C 1, we can approximate the quantity ^* 
in the expression for (see the second formula in ([57]) ) with 1. Thus ([57]) can be written as 

r* = f*^*, nt^ht + C- (38) 

Using (j27p . the posterior pdf for r* and rif is related to P{£,* , C*) iii tbe following way 

Pir*,nt) = ^P{e,n ^ ^P{C,C)- (39) 

Note that, in the current approximation, the pdf for variables r* and nt has a bivariate normal 
form. 

Based on the above pdf (j39p . we can now evaluate the maximum likelihood estimators, standard 
deviations, and the correlation coefficient for the variables r* and nt- For the maximum likelihood 
values we get 

ntML -nt + (,ML -nt+ J2^J2^(aYb*)2 ■ 

The standard deviation are given by 



Ant ^ AC = l/^JJ:iJ:Yi4K)^■ 
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(41) 



Finally, it can be shown that the correlation between r* and nt vanishes 

cov(r*, nt) 



0. 



Y^cov(r*, r*)cov(nt, rit] 
It is interesting to point out that these results are consistent with the results in 



(42) 



401]. The con- 



straints, presented in this paper, on the tensor-to-scalar ratio r are exactly the same as those in 



40]. Unlike the analysis in the current paper. 



401 ] works with a single free parameter r and does 



not consider nt as an independent free parameter. 



D. Constraints on parameter r 

In Subsection IIII Cl using the posterior probability function P{r*,nt), we have investigated 
the tensor-to-scalar ratio r* and the spectral index nt defined at the best pivot wavenumber k* 
(corresponding to the best pivot multipole it). We can now proceed to the analysis of the tensor-to- 
scalar ratio r, defined at an arbitrary pivot wavenumber ko, and determine the possible constraints 
on this parameter. 

From the Eq.Q, we can express the tensor-to-scalar ratio r, in terms of r*, kt and the spectral 
indices nt and n^, in the following form 

Inr = Inr* + ntln(A;o/A;;) + {1 - ns)ln{ko/k^). (43) 

It can be seen that, for a fixed value of the spectral index (see ([9])), r depends on the parameters 
r* and nt- Thus, the properties of r can be determined using the posterior pdf P{r*,nt), which 
was analyzed in detail in Section IIII Cl In the case ko ^ k^ it is more illustrative to consider the 
variable Inr instead of the variable r. For this reason, when dealing with the maximum likelihood 
estimators of tensor-to-scalar ratio defined at pivot scale different from the best-pivot scale, we 
shall use the corresponding logarithms 

InrML = Inrjfi + utML Mh/kt) + (1 " ^^s) ln{ko/kl), (44) 

where r|,j^ and ntML ^^e the maximum likelihood estimators expressible in terms of the input 
parameters r*, fit and the data {Df^ (see (j40p ). The uncertainty of r can be expressed in terms 
of the uncertainties Ar* and Arif determined in (j41|) . leading to the following expression 

Alnr ~ (Ar*/f*)'^ + {ln{ko/k^)Ant)'^ , 

= ^(4*)2 + (ln(/co/fc*)C)2. (45) 
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The quantities ^* and entering the above expression, can be expressed through CMB power 
spectra due to RGWs Cf^ using and ([28]) . 

From (|45|) it fohows that Ar/r > Ar*/r*, with the equahty holding for k^. Thus, the 

smahest uncertainty on tensor-to-scalar ratio r is achieved for the choice of the pivot scale ko = k*. 
This justifies the title "best" pivot wavenumber for k^. For a choice of pivot wavenumber ko ^ k^ 
the uncertainty in determining r becomes larger due to the uncertainty in determining the spectral 
index nt- 

Although, as was shown in Subsection IIII CI the quantities r* and nt are uncorrelated, this is 
not true for the quantities r and rit in general. In order to describe the correlation between r and 
nt, it is convenient to introduce the correlation coefficient 

_ cov(nt,lnr) 

P{nt,\nr) = , ^ n 1 ^ ' ^ ^ 

y'cov(nt, jcov(ln r, mrj 
where the notation cov(-, •) for the covariance was defined in (|35p. Using this definition, along with 

(|43p and ()42p . the terms entering the above expression can be evaluated as 

cov(lnr, Inr) = (Alnr)^, 
cov{nt,nt) = (Ant)2, 

cov(ni,lnr) = cov{nt,\nr*) + {\n{ko/k*))cov{nt,nt) = {\n{ko/k*)){Ant)'^ . 
Taking into account (j4ip . the correlation coefficient can be presented in the following form 



_ / cf{Hko/k*t)r 

As expected, for choice the k^ = k^, i.e. when the pivot wavenumber is chosen at the value of 
the best pivot wavenumber, the correlation between r and nt vanishes. On the other hand, for 
I \n(kQ/kf)\ ^ 1, i.e. for values of the pivot wavenumber significantly different from the best pivot 
wavenumber, the correlation coefficient approaches unity, implying a strong correlation between r 
and nt- 



E. Statistical properties of maximum likelihood estimators 

The exact numerical values of the maximum likelihood (ML) estimators CX/lj Cml' '^ml' ''^tML 
and InrML discussed in the previous subsections depend on the CMB data {DJ}. Since the 
set {DJ} is a single realization of an underlying random process characterized by the pdf ()Aip . 
the precise values of the maximum likelihood estimators will depend on this realization. For 
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this reason, it is instructive to analyze the distribution of these maximum hkehhood estimators in 
various reahzations of the underlying random process specified by the pdf for estimators of the CMB 
power spectrum {Df}. Heuristically speaking, the mean value of this distribution characterizes the 
typical value for the ML estimators that we are likely to observe (for a specific input cosmological 
model), while the standard deviation characterizes the typical departure from the mean value. 

Let us firstly, for simplicity, consider the estimators ^^^-^^ and Cml- '^^^ expectation values for 
these estimators can be calculated in the following manner 

(eW = ((E.Ey^r^D/lE.Ey aD) = {ZeEv 4 {4)) / {ZeZy 4') = h 

(Cml) = {{EiEYaJdJb})/{Z,EYiaJb}n = {nZy 4 {dj)b}) / {Z.EviajKr) = 0. 

The angle brackets (...), in the above expression and elsewhere in the text, denote the ensemble 
average over the joint pdf ()Aip . Furthermore, in this pdf, the input values for the tensor-to-scalar 
ratio and spectral index are chosen as r = r and rit = fit respectively. In deriving the above 
expressions we have firstly used and ([28l) . We have also used the identity {df) = af which 
follows directly from ([5]), (jl4p and (|20p . Finally, in the bottom line, we have used the definition of 
the best pivot multipole (f30ll . Similarly, the standard deviations can be calculated to yield 

^«M.=C, %*„=C. (49) 

Proceeding in an identical manner, the expectation values and standard deviations for the 
maximum likelihood estimators J'J/^, TitML ^i^i^ML are given by 

{r*ML)=r*{CML)=r*, 

{ntML) = fit + {Cml) = (^0) 
{Inr ml) =hif* + {fit -Us + l) ln{ko/kt) = Inr, 

and 



As expected, from expression ()50p it can be seen that the constructed ML estimators are unbiased. 
Furthermore, the standard deviation of the estimator cri^rML strongly depends on the choice of the 
pivot multipole ko, and is minimal for the choice ko = k^- We shall numerically verify these results 
in the following section. 
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F. The dependence of results on cosmological parameters and experimental noises 



Let us now address the question of detection of RGWs in various CMB experiments. In order 
to quantify the ability to detect the sig nature of RGWs in the CMB data, it is convenient to define 



the signal-to-noise ratio as follows 



21, 



4C| 



S/N^—. (52) 



Using expression (j4ip we arrive at an elegant expression for the signal-to-noise ratio 

S/N 



Thus, the signal-to-noise ratio contains contributions from individual power spectra and individual 
multipoles. These contributions have a clear physical meaning. For a particular power spectrum 
and a particular multipole, they represent the ratio of the expected signal due to RGWs to the 
overall uncertainty. 

As was mentioned in Section [III B[ for the analytical estimations, we had assumed ^* ^ ^* (see 
(|32p ). We can now relate this condition to the requirement on the value of the signal-to-noise ratio 
S/N. Using Eqs.dHD, dlQ]) and dH]), we find that 



S/N. (54) 



^* Ar* Ar* 

Hence, the condition ^* » ^* corresponds to the requirement S/N » 1, i.e. to the requirement 
that the RGW signal may be well determined at a high signal-to-noise ratio. 

In the discussion above we have mentioned that the best pivot multipole i^, the signal-to-noise 
ratio S/N and the uncertainty in determination of the RGW spectral index Ant depend on the 
input cosmological model and the specifics of the CMB experiment. Let us analyze this dependence 
in more detail. 

The input cosmological model is determined by specifying the background cosmological model, 
along with the parameters determining the density perturbations and gravitational waves. The 
background cosmological parameters and contribution from density perturbations are fairly well 
constrained by the current observations 



38l ]. The variation of these parameters within the margin 
allowed by these constraints will not significantly alter our results. For this reason, we shall fix the 
background cosmological model using the values of the typical ACDM model ([8]) . We shall also fix 
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FIG. 2: The figures show the value of the best pivot multipole (left panel), signal-to-noise ratio 
S / N (middle panel) and the uncertainty in the RGW spectral index Ant (right panel) as functions 
of the tensor-to-scalar ratio r*. The solid lines correspond to the Planck instrumental noises (see 
([15])); the dashed lines correspond to noises from cosmic lensing (see (fT6]) ); and the dotted lines 
correspond to reduced cosmic lensing noise (see p^ ). 

the contribution of density perturbations at a value ([9]). Furthermore, numerical calculations show 
that the dependence of various parameters on the input value of the spectral index ht is weak, 
for this reason in evaluations of this section we shall set fit = 0. Thus, we shall be interested in 
the dependence of the parameters on value of the input tensor-to-scalar ratio r. FIG. [2] shows the 
values of quantities £j , S/N and Ant as functions of f*, calculated using the expressions (130]) . (f53ll 
and (liTD . 

As was explained in Section [III the specifics of the CMB experiment are determined by the 
noise power spectra, the cut sky factor and window function. In this section we shall consider the 
parameters l^, S/N and Ant for the three cases specified in Section |TT] (see and ([IS]))- The 
different curves (solid, dashed and dotted) on the three panels in FIG. [2] show the corresponding 
values of quantities it, S/N and Ant for these three noise scenarios. 

The left panel of FIG. [2] shows the best pivot multipole It as a function of the input tensor-to- 
scalar ratio f* which is defined with respect to the best pivot multipole. It can be seen that, for 
small values of r*, the best pivot multipole it is small. This behaviour can be easily understood. 
For small values of f*, the constraints on r* and nt mainly come from i?-mode power spectrum 



20 



40|. However, in the 5-mode the main contribution to the signal comes from large angular scales 
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40|]. Thus, for 



corresponding to £ < 10, where the signal is mainly due to cosmic reionization 
small r* , the constraints on parameters r and rif are most stringently determined at large angular 
scales corresponding to multipoles i < 10. For this reason, for small values of r* the best pivot 
multipole is small, corresponding to the scale at which the parameters r and nt are most stringently 
determined. On the other hand, for large values f*, the best pivot multipole also becomes large. 
This happens due to two reasons. Firstly, with an increase in value of r*, the relative contribution 
of the reionization contribution to the to the S/N decreases, while the relative contribution of the 
multipoles around ^ ~ 90 (where the -B-mode spectrum is expected to have a maximum) increases 
(see FIG. ED. Thus, the contribution of RGWs at higher multipoles {i ~ 100) becomes significant, 
which in turn increases the value of the best pivot multipole. Secondly, when r* is large, the 
contributions from the C, T, E power spectra become important in constraining r and nt [40] . For 
these power spectra, the main contribution to the signal comes from the multipoles 10 < £ < 100 



40]. This again leads to an increase in the value of i^. 

The middle panel in FIG. [2]shows the signal-to-noise ratio S/N as a function of f*. As expected, 
the signal-to-noise ratio rises with the increase of f*. Setting the threshold value of S/N = 2, we 
can determine the detection possibilities for the three considered examples: r* > 0.05 for Planck 
noises; f* > 1.5 x 10^^ for the case with cosmic lensing; r* > 3.7 x 10^^ for the case with the 
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42 



reduced cosmic lensing. These estimations are consistent with previous results 

Finally, the right panel in FIG. [2] presents the achievable constraints on the spectral index Ant 
as a function of f*. As expected, the uncertainty in determining the spectral index drops with the 
increase of the input value f*. For the case of Planck mission, the uncertainty in estimation of nt 
always remains fairly large. Even for large value f* = 1 the constraint on the spectral index is 
Ant = 0.08 (for comparison, the Planck mission will be able to achieve constraint of Ang = 0.0045 
on the spectral index of density perturbations [l|]). For a value f* = 0.1, the constraint on the 
spectral index is Ant = 0.25, which is too large to constrain inflationary models or to verify the 
consistency relation. Potentially, in an idealized situation with reduced cosmic lensing, for r* = 0.1, 
we can constrain the spectral index to the level Ant = 0.007. If this accuracy can be achieved in 
the future, it will place a fairly tight constraint on inflationary models. 



21 



10 100 1000 

Multjpole: I 

FIG. 3: The comparison of and (T£,s (which enter the expression for signal-to-noise ratio S/N 
(|53|) ). for models with f* = 0.3 (black) and f* = 0.7 (red). The solid lines show the 'signal'- 
term (i.e. power spectrum) £{i + l)Cf /2t:{ixK?), and dotted lines show the 'noise'-term £{1 + 
l)a£,s/27r(/iK^). The quantity cfj^b was calculated using the Planck noises (fTSl) . 

IV. COMPARISON WITH NUMERICAL SIMULATIONS 

In Section mil we have analytically studied the likelihood analysis of the RGW parameters r*, nj 
and r, as well as introduced the best pivot multipole (corresponding to the best pivot wavenumber 
k^) and explained its significance. We have analytically derived expressions for the uncertainties 
of the RGW parameters and the value of the best pivot multipole, in terms of the CMB power 
spectra, experimental uncertainties and the estimators of the CMB power spectra. In this section 
we shall compare the analytical results of the previous section with numerical simulations. We 
shall show that, although we have used a number of approximations, the analytical results are in 
good agreement with the exact numerical results based on the analysis of simulated data. 

This section is separated into two parts. In the first subsection, using a single simulated data set 
{D^}, we shall use the Markov Chain Monte Carlo (MCMC) techniques to construct the posterior 
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pdf for the RGW parameters. We shall calculate the uncertainties and correlations associated 
with the parameters, and compare these values with the analytical predictions in Sections IIII CI 
and IIII Dl In the second subsection we shall generate 300 samples of data sets {D^}. For each 
individual sample, using the posterior pdf P{r*,nt) we shall calculate the estimates for the RGW 
parameters r'^^j^, ^tML ^ML- Analyzing the distribution of these estimates, we shall evaluate 
the mean values and the standard deviation, and compare these with the analytical predictions 
from Section [III E[ 



A. Likelihood analysis of a single simulated data set 

In this subsection, from a single simulated data set {Df}, using the likelihood analysis pro- 
cedure, we shall derive the constraints on the tensor-to-scalar ratio and the RGWs primordial 
spectral index. 

In order to simulate the CMB data, we shall randomly draw a data set {Df}, from an underlying 
pdf (|Aip (see Appendix [A|) . This pdf depends on the input cosmological model and characteristics 
of the CMB experiment. We shall choose as an input cosmological model, a model with the 
background cosmological parameters given in ([5]) and the contribution of density perturbation Q- 
The input parameters for the RGW field will be chosen as 

r = 0.3, ht = 0.0. (55) 

To characterize the properties of the CMB experiment, namely the power spectra of noises, the 
cut sky factor and the window function, we shall adopt the values specified for the Planck satellite 
mission ([15]) [l]. 

In order to simulate and analyze the data, we proceed as follows: 

1) We generate a single data sample {L'I'ly = C, T, E, B; i = 2,3, ■ ■ ■, 1000}, drawn from the pdf 

m- 

2) Using ([5n|) . we calculate the best pivot multipole scale £j = 21.1 (corresponding to the best 
pivot wavenumber = 0.002Mpc~^). Note that, the value of does not depend on the concrete 
realization generated in Step 1. 

3) Using the MCMC method (see 



45l | for details), we construct the likelihood function £ as a 
function of two free parameters r* and nt, with the other cosmological parameters fixed at their 
"best-fit" values given by ([8]) and Q. Choosing a uniform prior we build the posterior pdf P{r*,nt) 
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FIG. 4: 2-dimension and 1-dimension posterior constraints for parameters: r* and nt (upper 
panels), and for parameters: Inr and nt (lower panels). The blue '+' in the left panels indicate the 
value of the input model parameters. 

(which is exactly equal to the likelihood function £ (see (fM|) and ([25|) )). 

4) Using the posterior pdf P{r*,nt), we find the maximum likelihood values {f'p^j^,ntML)j plot 
the contours corresponding to 68.3% and 95.4% confidence interval regions in the {r*,nt) plane 
surrounding these values. We also calculate the 1-dimensional posterior pdfs for variables r* and 
nt- From P(r*,nt), we calculate the uncertainties Ar* and Ant- Using the importance sample 
technique (see [45], |46[), we evaluate the correlation coefficient P{r*,nt) defined in (P2]) . 

5) We now choose a different value of the pivot wavenumber /cq = 0.05Mpc~^, corresponding to 
the value for the pivot multipole = 500. Using (j43p . we calculate the tensor-to-scalar ratio for 
this pivot wavenumber r as a function of the parameters r* and nt. From the posterior probability 
function P{r*,nt), using the importance sample technique, we can obtain the uncertainty A Inr 
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and the correlation coefficient P(nt,inr) defined in (06]). 

The results of the simulation and analysis is shown in FIG. |H The panels on the top show the 
constraints in the r* — rit plane (top-left), and the 1-dimensional posterior pdfs for r* (top-middle) 
and rit (top-right). The constraint on the parameters r* and nt, together with the correlation 
coefficient are as follows 

r* = 0.343+°;°^^, (68.3%C.L.); nj = -0.067+°; (68.3%C.L.); = -0.02. (simulation result4^6) 

For comparison, the analytical formulae (j40p . (j4ip and ()42p yield the following results for these 
quantities 

r^^ ± Ar* = 0.345 ± 0.047; tuml ± Anj = -0.062 ± 0.135; p(r-',nt) = 0- (analytical result!^^7) 

As can be seen, the analytical results (j57p are in good agreement with results of simulation (j56p . 

The bottom panels in FIG. H] show the constraints in the Inr — plane (bottom- left), and 
the 1-dimensional posterior pdfs for Inr (bottom-middle). As expected, the confidence interval 
contours in the Inr — nj indicate a strong correlation between Inr and nt- The corresponding 
constraints and correlation coefficient follows 

Inr = -1.299+°;^^^, (68.3%C.L.); = 0.95. (simulation results) (58) 

The analytical expressions (HH), (05]) and (jl7]) . yield the following results for these quantities 

In rML ± Alnr = -1.307 ± 0.456; P(nt,inr) = 0-94. (analytical results) (59) 

Once again, we find that analytical and the exact results are consistent with each other. 

Furthermore, we have applied the same simulation and analysis procedure to the case with the 
"cosmic lensing" noises (see (|16p ). for input values of tensor-to-scalar ratio f* = 0.1, 0.2 and 0.3. 
We found that in all these cases, the numerical estimations for Ar*, Alnr and Ant agree with the 
analytical expression to within 20%. Thus the analytical formulae for Ar*, Alnr and Ant seem to 
be accurate. 

B. Maximum likelihood analysis in numerous data simulations 

In this subsection, we shall discuss the distribution of the maximum likelihood estimators for 
the RGW parameters r|,j^, vml and ntML in multiple realizations. We shall generate a simulated 
CMB data set {Df} a number of times. For each individual realization we shall calculate the 
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FIG. 5: The values of the ML estimators from 300 simulations are shown projected onto the 
'^tML ~ '^ML plane (left panel), and utML — ^^i"ml plane (right panel). The red '+' indicate the 
value of the input model parameters. 

estimators r]^fi, tml and ritML- We shall then analyze the distribution of these parameters and 
compare these results with analytical calculations. 

In order to generate and analyze the data, we proceed in the following manner: 

1) A collection of 300 samples of data sets {L'|'|y = T, E, B,C;£ = 2,3,- • -,1000} is randomly 
generated from an underlying pdf f{Df,Dj,Df,D^), given in ()Aip . The input cosmological 
model and the noise characteristics of the CMB experiment are chosen in the same manner as in 
Section HVAl 

2) Using (i30]l . we calculate the best pivot multipole £1 = 21.1 (corresponding to the best pivot 
wavenumber = 0.002Mpc"^). Note that, the value of does not depend on the concrete 
realization generated in Step 1. 

3) For each individual sample, we construct the likelihood function £ as a function of variables 
r* and rif, which is equal to the posterior pdf P{r*,nt) (see ()24p and (j25p )). For each individual 



sample, an automated search (which uses the numerical technique of the simulated annealing 



m 
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determines the maximum likelihood estimators r|,j^ and ntML i^^ which the posterior pdf P{r*,nt) 
reaches a maximum). The calculated values rj^^ and ntML &re plotted in FIG. [5] (left panel). 
4) We adopt a different pivot wavenumber kQ = 0.05Mpc~^, corresponding to the value for the 
pivot multipole Iq = 500. From the set of values {r^jj^intML}-, we calculate the corresponding 
values of tensor-to-scalar ratio for the new pivot wavenumber r^v/L using (j43p . The resulting values 
are illustrated in FIG. [5] (right panel). 

The mean values and standard deviations for the quantities r|^j^ and ritML (shown in FIG. O 
(left panel)), obtained from the analysis of the simulated data, are 

{rl^i) ± o-rj,^ = 0.298 ± 0.046; {ntML} ± (^ntML = -0.001 ± 0.152. (simulation results) (60) 

For comparison, we can calculate the corresponding quantities using the analytical expressions 
derived in Section lllli Using ([50|) and ([5T|) , we obtain 

{rli^) ± Uri^^ = 0.300 ± 0.047; {utML) ± cfntML = ± 0.135. (analytical results) (61) 

Comparing (j6ip with (j60p . we find that the analytical expressions are in good agreement with 
results of numerical simulation. 

In a similar fashion, for the mean values and the standard deviation of quantity tml (shown in 
FIG. [5] (right panel)), we obtain 

([nrML) icTinrA/L ~ —1-252 lb 0.508. (simulation result) (62) 
The analytical expressions (j50p and (j5ip yield the following results 

(Inrjv/L) ± '^inrji/i = -1-204 it 0.456. (analytical result) (63) 
Comparing (j62p with (j63p . we find a reasonable agreement to within 10%. 

V. CONCLUSION 

In this paper, we have analyzed the potential joint constraints on the two parameters character- 
izing the RGW background, the tensor-to-scalar ratio r and the tensor primordial spectral index 
nt, achievable by the upcoming CMB observations. We have shown that, in general, there exists a 
correlation between the parameters r and n^. However, when considering the tensor-to-scalar ratio 
r* defined at the best pivot multipole number £*, the correlation between r* and nt disappears. 
Furthermore, the uncertainty Ar* has the least possible value. We have derived analytical formulae 
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for calculating i*, Ar*, Anj, Ar, and the correlation coefficient between r and n^. Using numerical 
simulations of future CMB data we have verified the robustness of our analytical estimations and 
have shown that our fairly simple analytical expressions agree with exact numerical evaluations to 
within 20%. We have also discussed the dependence of our results on the background cosmological 
model, the amplitude of the RGWs, and the characteristics of the CMB experiment. We have 
studied the dependence of the signal-to-noise ratio S/N along with the value of the best pivot 
multipole and the uncertainty Ant on the amplitude of the RGWs. We show that, although the 
Planck satellite will potentially be able to measure the tensor-to-scalar ratio to a level r > 0.05 (at 
2a C.L.), the uncertainty in determining the spectral index will remain fairly large Ant ^ 0.25 (for 
r = 0.1). Thus, for example, the Planck satellite will not be able to verify the so-called consistency 
relation rij = — r/8. In an idealized scenario, where the noises are limited by reduced cosmic lensing 
noise, the precision Ant ^ 0.007 (for r = 0.1) is achievable, thus, potentially allowing tight con- 
straints on possible inflationary scenarios. The analytical results presented here provide a simple 
and quick method to investigate the ability of the future CMB observations to detect RGWs. 
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APPENDIX A: EXACT PROBABILITY DENSITY FUNCTIONS FOR DJ AND 
LIKELIHOOD FUNCTION 



In 



22( 1 (see also 



id, 



48 



49]) we have derived the pdfs for the best unbiased estimators DJ of 
the various CMB power spectra Cj . These were derived under the assumption that the primordial 
perturbations (density perturbations and RGWs) are isotropic and homogeneous Gaussian random 
fields, and that the noises associated with the CMB measurements can be assumed Gaussian. In 
this appendix we shall briefly list the main results that have been used in the present paper. 
The joint pdf for the estimators DJ, Df, Df and has the following form 



f{Df,Dj,Df,Df) = f{Df,Dj,Df)f{Df), 



(Al) 
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where the pdf /(-D?) has the form of the distribution, 



2"-<=/^r(ne/2)(fjf )^ 

and the joint pdf f{Df, DJ, Df) is the Wishart distribution 

^ {"2{T3^ ((^ + - ^) } • 



In the above expressions (lA2ll and (|A3l), are the corresponding CMB power spectra, Nf are 
the noise power spectra, and Wi is the window function. The quantity = {2£ + l)/sky is the 
effective degree of freedom for a particular multipole i in the case of partial sky coverage with the 
cut sky factor /sky The quantities v, x, y, z are defined as follows 



V ^ n,{DfWf + N^)l{CfWf + iVf ), 



X = n^iDjWf + Nj), y = n^{DfWf + iVf ), z = n^D^Wf. 

In (|A2p . af is the standard deviation for the multipole coefficient af^. The quantities aj , af 
and in (|A3P are correspondingly the standard deviations and the correlation coefficient for the 
multipole coefficients aj^ and af^. These are expressible in terms of the CMB and noise power 
spectra in the following form 



aj = xlcJWf + Nj, af = JcfWf + NP, af = JcfWf + N^, 



Finally, the likelihood function £ introduced in Section [III Al is. up to a constant of normalization, 
the product of the joint pdf /(ZJf , DJ , Df , Df), i.e. 

C^llf{Df,Dj,Df,Df). (A4) 
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